Optical quantum memory for polarization qubits with V-type three-level atoms 



D. Viscor, 1 A. Ferraro, 1 Yu. Loiko, 1,2 R. Corbalan, 1 J. Mompart, 1 and V. Ahufinger 1,3 

1 Departament de Fisica, Umversitat Autdnoma de Barcelona, E-08193 Bellaterra, Spam 
2 Institute of Physics, National Academy of Sciences of Belarus, Nezalezhnasty Ave. 68, 220072 Minsk, Belarus 
3 ICREA - Institucio Catalana de Recerca i Estudis Avancats, Barcelona, Spam 
(Dated: December 30, 2011) 

We investigate an optical quantum memory scheme with U-type three-level atoms based on the 
controlled reversible inhomogeneous broadening (CRIB) technique. We theoretically show the pos- 
sibility to store and retrieve a weak light pulse interacting with the two optical transitions of the 
system. This scheme implements a quantum memory for a polarization qubit — a single photon in 
an arbitrary polarization state — without the need of two spatially separated two-level media, thus 
offering the advantage of experimental compactness overcoming the limitations due to mismatching 
and unequal efficiencies that can arise in spatially separated memories. The effects of a relative 
phase change between the atomic levels, as well as of phase noise due to, for example, the presence 
of spurious electric and magnetic fields are analyzed. 

PACS numbers: 42.50.Gy,42.50.Md,42.50.Ex 



I. INTRODUCTION 

Optical quantum memories aim at the storage and re- 
trieval of quantum states of light by means of physical 
systems with long-living coherences. The realization of 
high efficient quantum memories is an ambitious goal to- 
wards which there is a considerable worldwide activity 
[THZ]. Optical quantum memories are essential for quan- 
tum information processing: they play a central role in 
quantum computation with linear optics [8] and, more 
generally, in quantum networks [5]. Specifically, optical 
quantum communication suffers from degradation of in- 
formation along transmission at large distances. To cope 
with this, advanced strategies have been proposed re- 
quiring the use of quantum repeaters of which quantum 
memories are the core ingredient [TU] . 

In quantum communication with photons, logical 
qubits can be encoded in several ways, for example via 
polarization, time-bin, path, phase, or photon-number 
encodings. Concerning polarization, the logical |0) and 
|1) quantum states correspond to two orthogonal polar- 
ization degrees of freedom, for example, left- and right- 
circular polarizations. It is on the storage and retrieval 
of polarization qubits that we will focus on this work. 
Up to now, polarization-qubit memories have been ex- 
perimentally realized mainly in cold atomic ensembles 
and atomic vapors |llj . In these experiments, polariza- 
tion encoding is first transformed into path encoding and 
then stored into two spatially distinguishable ensembles 
(actually, two different memories). Retrieval is eventu- 
ally performed by mixing again the output pulses, thus 
transforming back to the original polarization encoding. 

The above mentioned realizations are attractive, how- 
ever they suffer from several practical drawbacks due to 
the extra step of splitting spatially the input state into 
two beams. As a matter of fact, in order to preserve the 
original input quantum state, the two retrieved beams 
must be indistinguishable when mixed back again (apart 
from what concerns their polarization). In particular, 



they must occupy the same spatial and temporal mode, 
implying perfect matching at the mixing stage. In addi- 
tion, the efficiency of the two memories must be the same, 
otherwise the two polarization components of the origi- 
nal beam would be retrieved in an unbalanced fashion. 
Such necessary requirements are experimentally demand- 
ing. Clearly, similar problems would affect also solid- 
state based memories, if two different solid-state memo- 
ries were used to store each polarization component along 
the lines of [T2]. Recently, the storage of polarization 
qubits without spatially splitting the input photon has 
been addressed [13l [14] . In [13] the qubit is stored in a 
single atom inside a cavity by means of stimulated Ra- 
man adiabatic passage whereas in |14j a dense ensemble 
of four-level atoms in a tripod configuration is investi- 
gated for storage and retrieval of polarization qubits. 

We show here that a quantum memory for polarization 
qubits can be implemented also in a simple three-level 
system without the extra step of splitting spatially the 
input state. In particular we base our analysis in the pho- 
ton echo technique [15] , The latter has been developed 
especially for the storage of time-bin qubits, being partic- 
ularly suited for broad-band pulses, in contrast with the 
proposal in 13J. However, we will show in this work that 
photon echo quantum memories could also offer a com- 
pact resource to store and retrieve polarization qubits. 

Based on photon echo processes, two main approaches 
have been developed for devising quantum memories in 
solids: controlled reversible inhomogeneous broadening 
(CRIB) 16 -24 and atomic frequency comb (AFC) quan- 
tum memories [5] [B] I25j . Here, we will mainly focus on 
the former, even if adaptation to AFC can be envisaged 
with a similar approach. Quantum memories based on 
CRIB were first proposed in a seminal article by Moi- 
seev and Kroll [16J, where an optical pulse was stored 
and retrieved in a Doppler broadened two-level atomic 
vapor. In this technique, since the medium is inhomo- 
geneously broadened, the pulse is absorbed collectively 
by all the atoms in a coherent way, so each one evolves 



2 



with a different phase depending on its detuning. In this 
scheme, an auxiliary ir pulse is used, after the absorption, 
to increase the storage time by transferring the atomic 
coherences to a long-living third level. Afterwards, an- 
other counter-propagating 7r pulse brings the population 
back to the original two-level transition and imprints into 
the medium the phase matching conditions needed for 
the retrieval of the light pulse in the backward direc- 
tion. As an automatic consequence of the Doppler effect, 
the inhomogeneous broadening of the original two-level 
transition is now reversed for the backward propagat- 
ing retrieved pulse. Thus, a backward propagating time- 
reversed replica of the original pulse is retrieved. Further, 
the CRIB technique has been extended to solid-state ma- 
terials |17j . The key point is that in appropriate materials 
the relevant two-level transition can be artificially broad- 
ened, and the broadening reversed in a controlled way 
after the absorption and storage of the incident pulse. 
In later contributions [THl [20] it was shown that, if the 
phase-matching operation is not performed, the retrieved 
pulse propagates in the forward direction. In this case, 
the re-absorption by the medium itself limits the output 
efficiency to ~ 54%. An alternative to this method is 
the so-called longitudinal CRIB or gradient-echo mem- 
ory (GEM) [21], in which the artificial inhomogeneous 
broadening is not created at each spatial point, but it is 
created in such a way that the detuning changes linearly 
along the propagation direction of the light. Thus, each 
frequency component of the incident pulse is absorbed 
at a different position inside the medium, and once the 
detuning is reversed, this allow for a perfect retrieval ef- 
ficiency even in the forward direction for large optical 
depths. In this work we focus in the original proposal, 
called the transverse CRIB in contraposition to the lon- 
gitudinal case, although the results shown here could be 
easily extended to the GEM scheme. 

In most of the previous works concerning CRIB, usu- 
ally only two-level systems have been considered — apart 
from the mentioned auxiliary third level used only to 
perform the phase-matching operation and increase the 
storage time. Here we extend the CRIB approach to 
a medium composed of three-level atoms, paying special 
attention to the influence of the phase between the differ- 
ent atomic levels in the storage and retrieval efficiencies. 
We consider a collection of inhomogeneously broadened 
V-type three-level atoms (see figure [I]) , initially prepared 
in their ground state. The medium interacts with a weak 
pulse with two orthogonal polarization components, that 
is to say, left- and right-circular polarizations, each one 
coupled with one of the two optical transitions. With this 
scheme we will show that a V-type three-level medium 
can be used as a quantum memory for single photons in 
an arbitrary superposition of polarization states, without 
the need of spatially splitting the input pulse. 

The paper is organized as follows. In Section [TT] we 
will present the physical system under investigation and 
derive the optical-Bloch equations that govern its dynam- 
ical evolution. We will also introduce the figures of merit 



that will allow us to assess the quality of the quantum 
memory protocol, namely the efficiency and the fidelity 
of the retrieval process. In Section |III| we will report 
the analytical solutions for the evolution equations of the 
system at the stages of pulse storage and retrieval for 
both backward and forward modes. We will show that 
under the weak excitation regime, the system is equiv- 
alent to two independent two-level quantum memories. 
Next, in Section IV we will analyze the effects — on both 



the efficiency and the fidelity of the quantum memory — 
of a relative phase change between the different atomic 
levels as well as of phase noise during the storage time. 
In Section [V] we will verify the analytical predictions by 
the numerical integration of the full set of optical-Bloch 
equations. Finally, in Section \VJ\ we will summarize the 
results and conclude. 



II. BASIC EQUATIONS OF CRIB FOR V-TYPE 
THREE-LEVEL ATOMS 

The scenario we address here is sketched in figure [l] 
A single arbitrarily polarized weak pulse propagates in 
the forward direction (+2) and interacts with a medium 
composed of V-type three-level atoms, whose transitions 
have been subjected to artificial transverse inhomoge- 
neous broadening (figurejlja)), induced for example with 
a transverse electric field gradient. Note that in order to 
have a complete absorption of the field the inhomoge- 
neous broadening width must be larger than the spectral 
width of the input pulse. We assume the same nomi- 
nal frequencies for the two optical transitions of the V- 
scheme and that the left- and right- circular polarization 
components of the weak pulse interact with a particular 
optical transition (figure [TJb)). 

In this situation, the temporal evolution of a single 
atom can be described by the following density matrix 
equations: 

d 

-au(z,t) = ia 13 (z,t)il* 13 (z,t) + c.c. (la) 



dt 

d_ 

dt 
d_ 

dt 



d_ 
dt 



d_ 
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0-22(2, t) = \a 23 (z, t)tt 23 (z, t) + c.c. (lb) 

0-12(2, t) = iuj 12 cr 12 (z, t) + icr 13 (z, t)n* 23 (z, t) 

-ia 32 {z ) t)il 13 (z ) t) (lc) 

0-13(2, t) = iuJi 3 a 13 (z, t) + ia 12 (z, t)Q, 23 (z, t) 

-i [0-33(2, t) - au(z,t)] n 13 (z,t) (Id) 

0-23(2,*) = iw 2 30 23 (2,i) + \a 21 {z,t)VL 13 {z,t) 



-i [033(21, t) - <r 22 (z, t)] ^23(2, t) (le) 

where an denotes the population of level Ujj is 
the atomic coherence between levels \i) and \j), fly = 

(dijEij^j /h is the Rabi frequency with Eij being the 
slowly varying electric field amplitude of the light compo- 
nent coupled with transition \i) O \j), d{j is the electric 
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FIG. 1. (a) Illustration of the physical system under investiga- 
tion: a single pulse with central frequency ujq and two circular 
polarization components enters a medium with transverse in- 
homogeneous broadening, e.g., by means of a transverse elec- 
tric field gradient, (b) Scheme of the V-type three-level atoms 
that compose the medium. The left- (right-) circularly polar- 
ized component of the incident pulse couples to transition 
1 3} «-> 1 2} (| 3} <-» |1}), f2i3 and f223 denote the corresponding 
Rabi frequencies and A13 and A23 the detunings. 

dipole moment of the corresponding transition, h is the 
reduced Planck constant, and w, 3 = — u>j is the transi- 
tion frequency between levels \i) and \ j). Note that, since 
we consider a closed system, the atomic populations must 
satisfy an + a 2 2 + (J33 = 1. In addition, for simplicity, 
we have not included any incoherent decay term in equa- 
tions ([I]). This is well justified since we restrict ourselves 
to the case for which the lifetimes of the excited levels 
and of the optical coherences are much longer than the 
total time of the storing and retrieving processes. Note 
that this is a common assumption done in the literature. 
However, spontaneous emission is not the unique deco- 
herence mechanism that could affect the efficiency or the 
fidelity of the quantum memory implementation. In par- 
ticular, in Section |IV| we consider different phase noise 
mechanisms affecting the atomic transitions due to, for 
instance, the transfer of the stored excitation back and 
forth from the excited state to the metastable one, or an 
imperfect isolation of the atomic ensemble from external 
electric or magnetic fields. 

In order to solve analytically equations 0, we gen- 
eralize the treatment detailed in [19] to the case of a 
three-level medium. With this aim, we split the opera- 
tors associated with the two dipole transitions as well as 
the Rabi frequencies of the two field components in for- 
ward and backward modes, denoted by the superscripts 
"f" and "b", respectively: 

* W (*> *) = t)e^t-k aZ ) + ^ t)e ^ot+k oZ ) 

(2a) 

% p {z,t) = ^ p (M)e i( " ot - fcoz) + Ql p (z,t)e i ^ t+k ^ 

(2b) 

where ujq is the central frequency of the pulse and fco the 
corresponding wavenumber. From now on, we set p = 3 
and p — 1,2. For a weak input pulse the dynamical 
evolution of the level populations can be neglected, and 
therefore, since the initial population is assumed to be in 
the ground level |3), one can assume that 033 ~ 1 and 
an = 022 = 012 ~ during the whole process. With 



these assumptions, and decomposing the field into for- 
ward and backward modes, equations ([I]) can be written 
as follows: 

+ia p M^(z,t), (3) 

where we have made explicit the dependence of the opti- 
cal coherence on the detuning A = W13 — loq — UJ23 — too, 
and p = 3, /i, v = 1, 2 with p =/= v. 

The propagation of the forward and backward modes of 
light, in a reference frame moving with the pulses (z — > z 
and t — ¥ t =p z/c), is given by the following equations: 

-fi£ p f (z, t) = T h W / G w (A) o%(z, t, A)dA , 

OZ J— co 

(4) 

where — (+) refers to the backward (forward) mode, 
G^p (A) are the inhomogeneous atomic frequency distri- 
butions and rjpp = g 2 Ndf ip /hc with g 2 = uj /2e V be- 
ing the coupling constant, £0 the vacuum permittivity, 
N the number of atoms in the quantization volume V, 
and c is the speed of light in vacuum. It is worth noting 
that, although the treatment is performed in the semi- 
classical formalism, the linearity of the equations ensures 
the validity of this model also at the single photon level. 
Thus, in what follows, both the atomic coherences and 
the fields could be interpreted as classical amplitudes as 
well as quantum operators. 

From ([3]) and Q, it is direct to see that the time 
reversed (t —> —t) equations for the forward propagat- 
ing modes are identical to the backward ones under a 
sign change in the detunings and in the field amplitudes. 
This symmetry in the optical-Bloch equations is indeed 
the basis for the CRIB protocol. Thus, once the forward 
propagating input light pulse has been completely ab- 
sorbed, the reversing of the detunings (A — » —A) should 
allow to retrieve a time reversed copy of the input pulse. 
This reversing operation can be achieved, for instance, 
by changing the polarity of the device that creates the 
inhomogeneous broadening. At the same time, one has 
to apply a position dependent phase matching operation 
to transform the forward components of the atomic ex- 
citations into backward components, such that the re- 
trieved field propagates in the backward direction. This 
phase matching operation can be performed, for instance, 
by applying two counter-propagating ir pulses to trans- 
fer the atomic coherences back and forth to an auxiliary 
metastable state [T7J [TS], which also allows for longer 
storage times [26 . If the phase matching is not per- 
formed, the atoms will reemit in the forward direction, 
process that is described by the forward modes of equa- 
tions Q and Q. 

In this work we study both the backward and forward 
retrieval protocols, whose quality is assessed by defining 
the retrieval efficiency and fidelity as follows. For each 
of the two field components involved in the process, the 
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efficiency is defined as: 

/do, 



Eff 



b,f 
MP 



^ipi z — z out, a;) 



/dw nj»(« = o,w) 



(5) 



where the superscripts f and b refer to the forward and 
backward protocols, and i7" p (0,cj) and f2^(z out , w) de- 
note the temporal Fourier transforms of the input and 
output fields, respectively. For the backward and forward 
retrieval schemes in a medium with length L, the output 
fields are evaluated at the medium surfaces z out = and 
Zout = L, respectively. Equation ([5| addresses the ener- 
getic balance of the storage and retrieval processes, that 
is to say, the storage and retrieval of each component 
of the field independently from the other one. However, 
when considering the storage of quantum fields, it is of 
relevance to focus on the phase between the two compo- 
nents (polarizations) of the retrieved qubit by comparing 
it with the phase of the input qubit. Since we are in- 
terested in the case of initial pure states, the fidelity is 
defined as the overlap between the initial state \ipia) an d 
the retrieved state, which in general, could be pure or 
mixed: 



l^inl^'^in)! 



(6) 



where g h,i is the density matrix of the backward and for- 
ward retrieved qubits. Note that for those cases in which 
the quantum memory is perfectly efficient the retrieved 
state is pure. This is due to the unitarity of the dynamics 
and the fact that the medium and the fields end up in 
a factorized state. On the contrary, for a non-unit effi- 
ciency, for instance for finite optical depths, where part of 
the incident light is lost, the output state will be mixed. 



III. QUANTUM STORAGE IN U-TYPE 
THREE-LEVEL ATOMS 



In this Section we solve analytically equations ^ and 
Q, obtained under the weak input pulse approximation 
and considering that the atomic population is initially in 
the ground state. Under the weak field approximation, 
the population of the excited levels as well as the coher- 
ence between them remains negligible during the whole 
storage and retrieval process. In such case the evolution 
of the two components of the pulse are uncoupled from 
each other and equations ^ can be solved for each field 
component separately. In other words, the system turns 
out to be equivalent to two independent two-level media. 
Thus, the standard treatment of CRIB in two-level me- 
dia can be applied: each pulse component is absorbed by 
the transition it is coupled with, and the information is 
stored in the respective optical atomic coherence. After- 
wards, each pulse can be retrieved as in standard CRIB, 
either in the backward or in the forward direction. Thus, 
this protocol avoids the need of any active spatial sepa- 
ration of the pulses to store them. Therefore, the use of 



CRIB in U-type three level atoms does not suffer from 
detrimental effects due to the mismatch in spatially sep- 
arated memories. For later purposes, let us now present 
the solution for the field in all the stages of the CRIB 
protocol following the lines of [T§] . 



A. Absorption 

Within the weak input pulse approximation and, there- 
fore, neglecting the dynamics of the two-photon coher- 
ence <7i2, the equations for the two field components de- 
couple and equations ([3])-(|4]) can be analytically solved 
inserting the solution of equations Q into equations Q . 
Using the notation introduced in equation ^ and Fourier 
transforming the result, one obtains: 



d_ 

Or: 



where we have defined 

/OO pOO 
GW>(A)/ <■ <■ 
- oo JO 



ItJT " lAT drdA. 



(7) 



(8) 



To get analytical insight, we consider symmetric 
transitions, that is to say, equal electric dipolcs 



'■13 



'2:i 



= dj and inhomogeneous broadening dis- 
tributions (Gi 3 (A) = G 23 (A) = G (A)). Under these 
assumptions it follows that Hiq(u) = H-^aiuj) = H(oj) 
and ?7i3 = 7^3 = rj. Therefore the solutions of equations 
read: 



fC p (z, W )=lC p (0,u;)exp 



(9) 



B. Backward and Forward Retrieval 

We study now the propagation of the retrieved light 
pulse in backward and forward directions caused by the 
sign change of the detunings, after the absorption stage, 
at time t = 0. Note that for the backward propagat- 
ing case, a phase matching operation is also needed for 
the retrieval which corresponds to a sign change in the 
field amplitudes. Therefore, by changing A — > —A and 
fl^'p — > -fftp f for the backward and A — s- -A for the 
forward retrieval in equations ([3]), equations can 
be solved following the same method as in the absorption 
stage, using the boundary conditions f2 pp (0,w) = and 

f2 pp (L,cj) = for the forward and backward re-emission 
processes, respectively. By defining 

/ + OO /"OO 
G w (-A) / e^e-^drdA, 
-00 Jo 



(10a) 



J M » = / G w (-A) / e'^e-'^drdA 



(10b) 
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the equations for field components associated with each 
transition are: 



+ J^M tf?Jz,-u>) 



(11) 



with + (— ) for the backward (forward) configurations. 
The analytical solutions for the output fields read: 



(12) 



n^(z,o;) = -7 b ' t MnS.(0,-w), 

where we have assumed symmetric transitions and we 
have defined the backward and forward retrieval coeffi- 
cients as: 



7 b M = 



F(u)+H(-u] 



7 (w) = zrjJ{uS)e zr ' 2 x 



x sinhc I zrj 



F(oj)-H(-l 



(13a) 



(13b) 



with sinhc(x) = sinh(a;)/x denoting the hyperbolic sinus 
cardinal function. When the spectral bandwidth of the 
pulse is smaller than that of the inhomogeneously broad- 
ened transitions, the functions defined in Q and ( 10 1 can 
be approximated as F(lj) ~ H{—lo) ~ J(w)/2 ~ w/2. 
Therefore the absorption coefficient becomes a{oS) — 
rjJ(u>) ~ rjTT and 7 b,f (u>) simplify to: 



7 



(1 



-aL 



) and 



7 ~ aLe 



-aL/2 



(14a) 
(14b) 



with aL being the optical depth. Equation (12 1 shows 



that for the backward scheme both field components are 
recovered up to the factor 7 b , which approaches unity for 
a large enough medium length. This result confirms that 
ideally both polarization components can be stored and 
recovered perfectly in the backward direction, thus per- 
mitting to build up a quantum memory for a polarization 
qubit. On the other hand, in the forward retrieval case 
the quantum memory yields a forward retrieval with a 
maximum achievable efficiency around 54% as expected 
from [19] . 



IV. EFFECTS OF PHASE NOISE 

We consider now the effects of a specific variation in the 
phase of the atomic coherences. This is propaedeutic to 
the analysis of the role of the phase noise in the quantum 
memory, which will be the main goal of this section. Let 
us consider that, after the complete absorption of the 
field and before the reversing of the detuning, the phase 
difference between the atomic levels can be artificially 
engineered, for instance, by means of external electric 
or magnetic fields. The phase variation of the atomic 
coherences can be written, in the most general case, as: 



an i y e 



(15) 



with i, j — 1, 2, 3. This phase variation affects the quan- 
tum memory protocol yielding an additional phase to the 
initial conditions used in the retrieval stage of the pro- 
tocol. As a consequence, equations (111 transform as 
follows: 



dz >" 



b.r 



± 



F^ p (u;)Q h /Jz,u;) 



+ J„{w)&*»-*JVt%{z,-u,)\ 7^,(16) 



with + (— ) corresponding to the backward (forward) re- 
trieval protocol. The solutions of equations (fl6| read: 



7^)^(0,-^). (17) 



The above expression is identical to equations ( 12 ), apart 
from a phase factor. As a consequence, the retrieval ef- 
ficiency for each field is identical to the standard CRIB 
protocol. However the additional phase change intro- 
duced here affects the fidelity, defined in equation (|6j). 
Let us focus on the most interesting case corresponding 
to the backward retrieval, for which each component of 
the pulse can be recovered with unit efficiency for large 
optical depth media (i.e., 7 b (w) — > 1). Consider a generic 
polarization qubit at the input: 



c L \L)+c R \R), 



(18) 



where \L) and \R) denote a left- and right- circularly 
polarization states of a single photon, respectively. In 



the case of unit efficiency, the output backward field (17) 
would be: 



c L ^ L \L)+c R ^ R \R) 



(19) 



being (f>L = 4>3 ~ 01 an d §r = 4>3 ~ 4>2- From equation 
(|6j), the corresponding input-output fidelity reads 



F 



ICR 



(20) 



with (p = cf>i — 4>2- Clearly, even if each component of 
the initial qubit is completely recovered, the additional 
phase accumulated in levels 1 and 2 degrades the fidelity. 
The latter is instead unaffected by a phase change in 
the lower level since, in this case, no phase difference is 
accumulated between the coherences ai3 and 023. 

The analysis performed up to here has been referred 
to a system evolving in the absence of noise. Let us con- 
sider now the presence of noise in the phase of the atomic 
transitions. Two mechanisms can be foreseen that yield 
phase noise in the atomic transitions during the storage 
time. First, consider that after the pulse is absorbed, 
each atom — apart from evolving according to its free 
Hamiltonian — is affected by spurious magnetic and elec- 
tric fields. This would lead to a non-ideal rephasing after 
the reversing of the detuning. Second, as it has been 
discussed in the Introduction, in order to increase the 
storage time it is customary to temporally transfer the 
population from the excited unstable level to an auxiliary 
level endowed with larger coherence time. This procedure 
is performed via the coupling with light pulses of area 7r 
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FIG. 2. (a) von Mises distribution of the phase of the lower 
level 1 3), for different values of the inverse width distribution: 
&3 = (solid line), ks = 1 (dashed line), = 5 (dot-dashed 
line), k% = 20 (dotted line). Memory efficiency for the (b) 
backward and (c) forward retrieval schemes as a function of 
the optical depth for different amounts of phase noise affect- 
ing the ground level after the absorption of the pulse. The 
different lines correspond to the same inverse widths of the 
distribution as in (a). 



[16HT8] . Ideally, these pulses should couple to the system 
without introducing any uncontrolled phase shift to the 
original atomic transition. However, in any realistic sit- 
uation, the intensity and phase fluctuations of the light 
pulses will result in a phase shift, probably different for 
each atom, when the population is transferred back to 
the original excited level. Both aforementioned noise 
mechanisms can be modeled by adding a phase noise to 
each atomic transition. In practice each atomic coher- 
ence is subjected to phase transformations as those given 
by equation (15). Since the additional phase is random 
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FIG. 3. Maximum memory efficiency for the backward (solid 
line) and the forward (dashed line) retrieval schemes as a 
function of the inverse width of the phase noise distribution 
k 3 . 



and centered in zero, we consider a von Mises phase dis- 
tribution (circular normal distribution), that is to say, a 
continuous probability distribution on a circle: 



Mi 



i) 



^kj cos <f>j 

27r/ (%; 



(21) 



In the definition above j refers to the level affected by the 
noise, <pj is the phase angle, Ia{kj) denotes the modified 
Bessel function of zero order, and kj is the inverse width 
of the phase noise distribution, that is assumed to be 
centered at zero phase value. For kj = the distribution 
is uniform along the interval [— tt, ir] and the larger the kj, 
the narrower the noise distribution, and vice versa (see 
figure ^ a)). Notice that the actual form of the chosen 
distribution is however not relevant for our purposes. 

In order to introduce the effects of phase noise in equa- 
tions ([3j)-(|4p, one can proceed as in equation (15), but 
now an additional averaging over all the possible phases 
must be performed in the field equations: 



d_ 
dz 



/7T />0O 
-7T J — OO 

x G MP (A) ( 7^ f (z,t,A> 1 > 2 ,0 3 )] • 



(22) 

where Mt = M\ (<j>i) (fa) M3 (^3). For a generic 
noise acting upon all the atomic levels, the expressions 
for the retrieved fields, both in the forward and back- 



ward configurations, read as in equations (17) except for 
a multiplying factor: 

f^ p f (z, W ) = ~K p K^ b > f (uj)^ p (0,-uj), (23) 

where - e ±i ^M M (^)d^ = /i(fc M )// (fc M ) that 
increases, monotonously with fe„, from zero to one. 

Consider, for instance, that the phase noise acts upon 
level 1 3) only (the case of noise acting on the other lev- 
els being similar) . From expressions ( 23 1 , the retrieval 
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efficiencies for the backward and forward modes can be 
obtained: 



Eff 



b.f 
MP 



(^37 



b,f\ 



(24) 



where 7 b and 7 f are defined in ( 14a| and (f4b), respec- 
tively. Notice that the efficiency does not depend on 
the amplitude components of the initial field, so it is in- 
dependent of the particular superposition state we are 
considering. Focusing in the backward scheme, we plot 
in figure [2jb) the retrieval efficiency as a function of the 
optical depth for different values of the inverse width of 
the distribution fc 3 . We see that for large fc 3 (narrow 
phase noise distribution) the efficiency is unaffected by 
the presence of noise, and a perfect recovery of the in- 
put light pulse can be reached for a large enough optical 
depth. However, as fc 3 becomes smaller (a wider noise 
distribution), the efficiency suffers from relevant detri- 
mental effects. In this case, the input pulse cannot be 
recovered completely, even for a large optical depth. For 
the forward retrieval scheme, similar results are obtained 
(see figure |2jc)) with the obvious difference that in this 
case the recovery efficiency is limited to ~ 54%. These 
results do not contradict the above mentioned fact that 
both the efficiency and fidelity are unaffected by a single 
specific phase change 3 of the ground state (see equa- 
tions (17), and ( 20 1 ) . In fact, when the phase noise is 
present, the situation is different: since each atom ac- 
cumulates a different phase during the storage time, the 
rephasing process, which is a collective phenomenon, can 
not be achieved properly. 

In order to have a compact assessment of the effect of 
the phase noise, we plot in figure [3] the maximum value 
of the backward (solid line) and forward (dashed line) 
efficiency as a function of the noise parameter fc 3 , for the 
case where the phase noise is only included in level |3). 
We can see that the memory efficiency in both backward 
and forward schemes degrades substantially as the phase 
noise distribution becomes wider. Note that the actual 
values of the noise depend on the chosen model, however 
the qualitative behavior is a general feature. 

In order to obtain a complete picture of the retrieval 
process in the presence of noise, the input-output fidelity 
has to be considered too. As indicated, whenever noise is 
present the retrieval efficiency is reduced with respect to 
the ideal case, being below the unity also in the backward 
scheme. This implies that the output state q°' is not 
pure any more. The expression of g°' can be obtained, 
however, simply by modeling the non-ideal memory with 
an ideal one preceded by a beam-splitter with transmit- 
tivity equal to the efficiency of the memory. The reflected 
beam is then traced out. In this way, one obtains that 
the fidelity is equal to the retrieval efficiency, in the case 
of polarization qubits. 



V. NUMERICAL RESULTS 

Analytical results reported in previous sections are 
based on the weak field approximation in which the tem- 
poral dynamics of the atomic populations and the two- 
photon coherence uyi was neglected. In this section we 
verify the validity of the analytical approach by perform- 
ing numerical integration of the full optical-Bloch equa- 
tions ([3]|4| with a finite difference method beyond this 
approximation. Atoms with closely spaced frequencies 
are grouped in equidistantly spaced frequency classes (la- 
beled by an integer number n) with central frequencies 



<jj 



(n) 
13 



J 23 



At each spatial point z the inhomogeneous 
spectral distribution (provided by the transverse broad- 
ening) of the frequency classes is given by the function 
G (A&A which, for simplicity, is assumed to be a rectan- 
gular distribution of width much larger than the spectral 
width of the input pulse. 

The calculation scheme is performed in the follow- 
ing steps: (i) From the temporal profile of the incident 
field pulse fi MP (z,i), we calculate the temporal evolu- 
tion of the density matrix at each spatial point z and 
for each frequency class by integration of equations ^ 
with time step At. (ii) This temporal evolution given 
by c^p (z,t, A^) is used in equations Q to propagate 
the field pulses between the spatial points z and z + Az. 
Further, steps (i) and (ii) are iterated. Temporal At 
and spatial Az steps are chosen small enough and the 
number of frequency classes large enough to assure the 
convergence of the calculations. 

During the absorption stage, the pulse propagation is 
initiated at z = 0, where the density matrices of all atom 
frequency classes are known. As in the analytical ap- 
proach, we assume that initially all atoms are in their 
ground level, i.e., cr 33 (z, t — — 15/xs) = I. The input field 
pulses at z = are assumed to have Gaussian temporal 
profiles and we define the normalized intensities of each 
component as 



fii 3 (o, t c ) r+|o 23 (o, t c ) r 



(25) 



where fi = 1,2, p = 3 and t c is the time corresponding 
to the peak of the pulse. The spatio-temporal dynam- 
ics of Ii 3 (z,t) and I 2 -i{z,t) during the absorption stage 
of the protocol is shown in figures ga) and (d), respec- 
tively. The initial intensities of the components are taken 
as h 3 (0,t c ) = 0.6 and I 2 3(0,t c ) = 0.4. Note that the ab- 
sence of field at the output surface, corresponding to an 
optical depth az — aL = 4.5 with L being the length 
of the medium, evidences the efficient absorption of both 
field components. 

For the simulation of the forward and backward re- 
trieval protocols we reverse the sign of the detuning for 
each frequency class A' n ) — > —A^ n '. For the backward 
protocol, we additionally change the sign of the field to 
implement the phase matching operation. After this step, 
the phase noise is introduced by imposing at each spatial 



Absorption Backward Forward 

(a) (b) (c) 
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FIG. 4. Contour plots of the spatio-temporal dynamics of the 
normalized intensities of the two circular polarization com- 
ponents of the field I\z(z,t) (a-c) and l23(z,t) (d-f) in an 
inhomogeneously broadened V-system with phase noise dis- 
tribution of width &3 = 5 in the lower level |3). Figures (a,d) 
show the absorption stage, while (b,e), and (c,f) correspond 
to the backward and forward retrieval schemes, respectively. 



position a phase distribution in the atomic coherences 
of each frequency class, n, given by equation (21), i.e., 

0) 



(n) [z,t 



0) 



"° . The ensemble av- 



;(«) _ / a (n) 



erage afcu — (crfcv J is used as the initial density matrix 

element for each frequency class at the retrieval stage. 
During the retrieval stage the field propagation is evalu- 
ated by means of equations ^ and Q from the surface 
az = {az = aL) to the surface az = aL (az — 0) for 
the forward (backward) protocol. The spatio-temporal 
dynamics of the retrieved field intensities, /13 and J23, 
in the presence of a phase noise distribution in the lower 
level 1 3) of inverse width k% — 5 are shown in figures |4](b) 
and (e) for the backward retrieval and in figures |4|c) and 
(f) for the forward retrieval protocols. 

For the parameters used in our simulations we have nu- 
merically observed that the two-photon coherence (C12), 
the population of the excited levels (<tii, 022), and the 
phase difference between the field components did not 
exceed 10 -6 during the whole dynamics, validating the 
weak field approximation used to derive the analytical 
solutions in Section |III| By comparing the absorption of 
the field components (figures Qa) and (d)) with the cor- 
responding retrieved fields (figures |4|b) , (c) and (e),(f), 
respectively), we note that, although the efficiency of the 
re-emission is lower than in the ideal case where no phase 
noise is present, the decreasing on the retrieval intensities 
is the same for both components. 

We have checked that beyond this particular case, the 
propagation dynamics and the maximal retrieval efficien- 
cies acquired from the numerical simulations for both 



forward and backward protocols are in good agreement 
with the analytical results, confirming the validity of 
the approximations that we assumed in the analytical 
approach. In particular, for the parameter values of 
figure [4j the maximum backward and forward numer- 



ically evaluated efficiencies read Eff 
Eff 



(num) 



= 0.79 and 



(num] 



= 0.43, respectively that match the analytical 



results shown in figure [2] for k% = 5 (dot-dashed line) . 



VI. CONCLUSIONS 

In this paper we have shown that optical quantum 
memories based on I^-type three-level atoms can be used 
to store and retrieve polarization qubits by means of 
the CRIB technique without the need of two separate 
quantum memories or spatially splitting the input state, 
hence presenting experimental compactness. In the weak 
field regime, the optical-Bloch equations for the three- 
level atoms decouple into two sets of equations describ- 
ing the evolution of the two circular polarization compo- 
nents of the field resulting in two independent two-level 
quantum memories. In this scenario, we have analyti- 
cally studied the effect of an extra phase in the atomic 
levels. We have shown that when the phase variations 
affect equally all the atoms of the sample, the fidelity 
is degraded whereas the efficiency for each polarization 
component of the pulse is unaffected. In contrast, when 
we consider an inhomogeneously distributed phase noise 
for the different levels, the rephasing process during the 
retrieval stage is not properly achieved and both the ef- 
ficiency and the fidelity are degraded. In addition, we 
have numerically checked the assumptions performed in 
the analytical approach, by numerically integrating the 
full optical-Bloch equations. The good agreement be- 
tween the numerical and analytical results demonstrates 
the validity of the analytical results obtained under the 
weak field approximation. 

Along this work we have proposed the implementation 
of a quantum memory for qubits in a superposition state 
of left and right circular polarizations in y-type three- 
level atoms. In the original CRIB proposal [18 rare- 
earth-ion-doped (REID) crystals were suggested as pos- 
sible candidates to implement quantum memories due to 
its well-known energy level structure and long coherence 
times at low temperatures [55]. However, although V- 
type three-level configurations have been implemented in 
REID crystals (see for instance [27] for a EIT experiment 
in Pr 3+ : Y 2 Si05), the optical transitions in this kind of 
systems are not polarization selective, so the only way to 
address them separately is to use a source whose band- 
width is smaller than the ground-state sublevel splitting. 
This restriction is an important drawback for the imple- 
mentation of our proposal in REID, where a system with 
specific selection rules is needed to absorb the left and 
right circular polarization components of the photonic 
qubit. Nevertheless, atomic vapours do not suffer from 
these limitations and could be good candidates to implc- 
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ment a polarization qubit three-level memory. In par- 
ticular, regarding photon echo quantum memories, warm 
[21] and cold [52] atoms of 87 Rb have been shown to be 
useful for A-type GEM implementations, with reasonable 
high efficiencies and storage times. 

Finally, we are presently investigating the implemen- 
tation of optical quantum memories based on CRIB for 
A-type three-level atoms. In this last case, the atomic 
coherence between the two ground states plays an im- 
portant role in both the absorption and retrieval stages 
and, in general, the system does not resemble two inde- 
pendent two-level quantum memories. The dynamics of 
this system, being much more involved and opening new 
scenarios for the manipulation of polarization qubits, is 



out of the scope of the present paper and will be discussed 
elsewhere. 
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